  use "$data\women_analysis.dta", clear

  global controlwoman mbase_ageyears  fdstotal rantime ranerror numberkids
  global controlhh land members adults mbase_electrified  mbase_OB mbase_scheduled

  estimates clear
  ****Table 6
  foreach i in score factualscore cognitivescore decision_making_index cook2 marketin2 marketout2 health2 education12 ///
  education22 status_index usevehicle personelmobile personelbankaccount jewellerywithyourself veil{
  ivregress 2sls `i' (participatedTA_IV=ITTanalysis) i.strata $controlwoman $controlhh if sequence==1, vce(robust)
  est store `i'
  su `i' if sequence==1 & ITTanalysis==0
  estadd scalar controlmean=`r(mean)'
  }
  
 esttab * using "$results\Table6.xls", stat(controlmean N r2) keep (participatedTA_IV) b(3) se(3) replace
 
 ***FDR q-values
 *The FDR Q-values are computed following  Benjamini and Hochberg (1995) within each category.
 
 exit
 






